LARGE DEVIATIONS FOR LOCAL TIME FRACTIONAL BROWNIAN 

MOTION AND APPLICATIONS 



Abstract. Let W'^ = {W'^{t),t e R} be a fractional Brownian motion of Hurst index 
H £ (0, 1) with values in R, and let L — {Lt,t > 0} be the local time process at zero of 
a strictly stable Levy process X = {Xt,t > 0} of index 1 < a < 2 independent of W'^ . 
The Qf-stable local time fractional Brownian motion — {Z^(t),t > 0} is defined by 
(t) = W^{Lt). The process Z^ is self-similar with self-similarity index H{1 — —) and 
is related to the scaling limit of a continuous time random walk with heavy-tailed waiting 
times between jumps ([SlISl]). However, Z^ does not have stationary increments and is 
non-Gaussian. 

In this paper we establish large deviation results for the process Z^ . As applications 
we derive upper bounds for the uniform modulus of continuity and the laws of the iterated 
logarithm for Z^ . 



Self-similar processes arise naturally in limit theorems of random walks and other sto- 
chastic processes, and they have been applied to model various phenomena in a wide range 
of scientific areas including telecommunications, turbulence, image processing and finance. 
The most important example of self-similar processes is fractional Brownian motion (fBm) 
which is a centered Gaussian process = {W"{t),t G M} with W"{0) = and covari- 
ance function 



where H G (0, 1) is a constant. By using (jl.ip one can verify that is self-similar with 
index H (i.e., for all constants c > 0, the processes {W^{ct),t E M} and {c^W^it),te M} 
have the same finite-dimensional distributions) and has stationary increments. When H = 
1/2, is a two-sided Brownian motion, which will be written as W. 

Many authors have constructed and investigated various classes of non-Gaussian self- 
similar processes. See, for example, |39j for information on self-similar stable processes 
with stationary increments. Burdzy [10^ ITT] introduced the so-called iterated Brownian 
motion (IBM) by replacing the time parameter in W by an independent one-dimensional 
Brownian motion B = {Bt,t > 0}. His work inspired many researchers to explore the 
connections between IBM (or other iterated processes) and PDEs [H [H HI [35] , to establish 
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potential theoretical results [SJ [THJ [331 EH 133 EZ] and to study its sample path properties 

In this paper, we consider another class of iterated self-similar processes which is related 
to continuous-time random walks considered in [GlEl]. Let = {W^{t),t G M} be a 
fractional Brownian motion of Hurst index H E (0, 1) with values in R. Let X = {Xt,t > 0} 
be a real-valued, strictly stable Levy process of index 1 < a < 2. We assume that X is 
independent of . Let L = {Lt, i > 0} be the local time process at zero of X (see Section 
2 for its definition). Let = {Z^{t),t > 0} be a real- valued stochastic process defined 
by Z^{t) = W^{Lt) for all t > 0. This iterated process will be called an a-stable local time 
H -fractional Brownian motion or simply local time fractional Brownian motion. 

Since the sample functions of and L are a.s. continuous, the local time fractional 
Brownian motion Z^ also has continuous sample functions. Moreover, by using the facts 
that is self-similar with index H and L is self-similar with index 1 — 1/a, one can 
readily verify that Z^ is self-similar with index H{\ — 1/a). However, Z^ is non-Gaussian, 
non-Markovian and does not have stationary increments. When H = 1/2, we will call Z^/"^ 
the local time Brownian motion and denote it by Z for convenience. 

The local time Brownian motion Z emerges as the scaling limit of a continuous time 
random walk with heavy-tailed waiting times between jumps [6l[3T]. Moreover, local time 
Brownian motion has a close connection to fractional partial differential equations. Baeumer 
and Meerschaert [3j showed that the process Z can be applied to provide a solution to the 
fractional Cauchy problem. More precisely, they proved that, if Lt is the local time at of a 
symmetric stable Levy process, then u{t,x) = E,x[f {W (Lt))] solves the following fractional 
in time PDE 

(1.2) ■^u{t,x) = Axu{t,x); u{0,x)=f{x), 

where /? = 1 — 1/a and d^g{t)/dt^ is the Caputo fractional derivative in time, which can be 
defined as the inverse Laplace transform of s^g{s) — s^~^g{0), where g{s) = e~'^^g{t)dt 
is the usual Laplace transform. Recently Baeumer, Meerschaert and Nane [4j further es- 
tablished the equivalence of the governing PDEs of W{Lt) and VF(|-Bt|) when a = 2 and 
/3 = 1/2. Here B = {Bt,t > 0} is another Brownian motion independent of W and X. 
The process Z has also appeared in the works of Borodin [F, "5], Ikeda and Watanabe [26j, 
Kasahara [27], and Papanicolaou et al. ^38j. In [17j, Csaki, Foldes and Revesz studied the 
Strassen type law of the iterated logarithm of Z{t) = W{Lt) when Lt is the local time at 
zero of a symmetric stable Levy process. 

For all H S (0,1) and a E (li2], a-stable local time ff-fractional Brownian motions 
form a new class of self-similar processes. It is natural to expect that they arise as scaling 
limit of continuous-time correlated random walks with heavy-tailed waiting times and, as 
such, they are potentially useful as stochastic models. Hence it is of interest in both theory 
and applications to investigate their probabilistic and analytic properties. Due to the non- 
Gaussian and non-Markovian nature of local time fractional Brownian motions, the existing 
theories on Markov and/or Gaussian processes can not be applied to them directly and some 
new tools will have to be developed. The literature on iterated Brownian motion mentioned 
above provides an instructive guideline for studying local time fractional Brownian motions. 
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The objective of the present paper is to estabUsh large deviation results for the local 
time fractional Brownian motion and apply them to study regularity properties of the 
sample paths of . We will consider the interesting problem of determining the domain 
of attraction of Z^ in a subsequent paper. 

The following Theorems 1 1 . 1 1 and 1 1 . 2 1 are our main results. 

Theorem 1.1. Lei Z^ = {Z^{t),t > 0} be an a-stable local time H -fractional Brownian 
motion with values in M and 2H < a. Then for every Borel set I? C 

(1.3) limsupt" «-2-ff logPU" --^H Z^{t) eD \ <- \n^A\{x) 
and 

(1.4) liminft" — 2ii logPU" — z"{t) eD \ >- inf A];(x), 

i^oo I J x<=D° 

where D and D° denote respectively the closure and interior of D and 

a-2H 

a + 2H f a-2H 



(1.5) ^^(-) = ^^i^^J 

In the above, Bi = Bi{H,a,x,t^) is the positive constant defined by 

2H 

a-2H HA1 



(1.6) Si 



^1 



2a l(l-i)"-i; 

\ / 

where Ai is the constant given by 

r(l - i)r(i)x^/"cos(iarctan(i/tan(^))) 



(1.7) 



7ra[l + (i/tan(^))2]i/(2") 



and u ^ [—1,1] and X > are the parameters of the stable Levy process X defined in \2. 

In (jl.5p and the sequel, for any 7 > and x G M, the term x^'^ is defined as {x^)^ . Since 
2H < a, one can see that the function A* in (jl.5p is even, convex and differentiable on M. 

2H(a-l) 2H(a-l) 

In the terminology of [20] , Theorem 11.11 states that the pair 

(t- c,-2H Z" {t),t «-2ff ) 

satisfies a large deviation principle with good rate function AJ. When H = 1/2, it yields a 
large deviation result for the local time Brownian motion and, moreover, the constants Bi 
in (|1.6p can be simplified. 

Letting D = [3;,oo), we derive from Theorem 11.11 and the self-similarity of Z^ the as- 
ymptotic tail probability P|Z^(1) > x} as x ^ 00. The following theorem is more general 
because it holds for all H G (0, 1) and a £ (1, 2]. 

Theorem 1.2. Let Z" = {Z^{t),t > 0} be an a-stable local time H -fractional Brownian 
motion with values in M. Then for any < a < b < 00, 

logF{\Z^{b)-Z"{a)\>x} 
(1.8) hm ' ' = -B2, 
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where B2 = B2{H,a,x,'^) is the positive constant defined by 



(1.9) .««. 

In order to prove Theorems 11.11 and II. 2[ we first study the analytic properties of the 
moment generating functions of Z^{t) and \Z^{b) — Z^{a)\. This is done by calculating 
the moments of (t) and \Z^ (b) — Z^{a)\ for < a < 6 directly and by using a theorem of 
Valiron |^40j. Then Theorems ll.ll and ll.2l follow respectively from the Gartner-Ellis Theorem 
(cf. [20] ) and a result of Davies [18] . 

The rest of the paper is organized as follows. In §2, we derive sharp estimates on the 
moments of the local time Lt of X, and the moments of Z^{b) — Z^{a). These esti- 
mates are applied in §3 to study the analyticity of the moment generating functions of 
Z^ , and to derive large time behavior of the logarithmic moment generating functions 
logE[exp((9Z^(t))] and logE[exp(t|Z^(6) - Z^(a)|^)] for suitably chosen /3 > 0. In § 4, 
we prove Theorem ll.il and II. 21 We will also establish similar tail estimates for the maximum 
maxjgjd \Z^{t) — Z^{a)\. In §5, by combining the large deviation result with the methods 
in [13\ Theorem 3.1], we establish local and uniform moduli of continuity for Z^ . We also 
obtain an upper bound in the law of the iterated logarithm for Z^ . 

Acknowledgment This paper is finished while Y. Xiao is visiting the Statistical & Applied 
Mathematical Sciences Institute (SAMSI). He thanks the staff of SAMSI for their support 
and the good working conditions. 

2. Moment estimates 

A Levy process X = {Xt,t > 0} with values in M is called strictly stable of index a G (0, 2] 
if its characteristic function is given by 

r,Y^^ ( .i,,. ! +^^sgn(g)tan(^) ^ 

(2.1) E[exp«At)J = exp I ^ 1 , 

where — 1 < < 1 and x > are constants. In the terminology of [39^ Definition 1.1.6], v 
and x"^''" respectively the skewness and scale parameters of the stable random variable 
Xi. When a = 2 and x = 2, X is Brownian motion. In general, many properties of stable 
Levy processes can be characterized by the parameters a, v and x- For a systematic account 
on Levy processes we refer to [7]. 

For any Borel set ICR, the occupation measure of A on / is defined by 

(2.2) ^ii{A) = \i{t £ I : Xt^A] 

for all Borel sets ACM, where Ai is the one-dimensional Lebesgue measure. If is 
absolutely continuous with respect to the Lebesgue measure Ai on M, we say that X has a 
local time on I and define its local time L{x,I) to be the Radon-Nikodym derivative of iij 
with respect to Ai, i.e., 

L{x,I) = ^{x), VxGM. 
uAi 
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In the above, x is the so-called space variable, and / is the time variable of the local time. 
If / = [0,t], we will write L(x,I) as L{x,t). Moreover, if x = then we will simply write 
L{0,t) as Lf 

By using a monotone class argument, one can verify that L{x,I) satisfies the following 
occupation density formula: For every measurable function / : R — > M_|_, 

(2.3) / f{X{t))dt= [ f{x)L{x,I)dx. 

J I JR 

It is well-known (see, e.g. [.7J) that a strictly stable Levy process X has a local time if 
and only if a G (1,2]. In the later case, L{x,t) has a version that is continuous in {x,t). 
Throughout this paper, we tacitly work with such a version so that the local time process 
L = {Lt,t > 0} has continuous sample paths. 

It follows from (12. and the self-similarity of X that L{x,t) has the following scaling 
property: For every constant c > 0, c^~^/"L(c~^/"x, c~^t) is a version of L(x,t). In partic- 
ular, by letting x = we see that Lt is self-similar with index 1 — — . 

For the purpose of the present paper, it will be convenient to express the local time L{x, t) 
as the inverse Fourier transform of ^{u,t) := fL^t]{u), namely 

L{x,t) = — / exp{—iux)Ji{u,t) du 

(2.4) -'^ 



— [ [ exp ( — iux + iuX{s)) duds. 
2vr Jo Jr 



This formal expression can be justified rigorously (see [21]). Moreover, it follows from (25.2) 
and (25.7) in [21j that for all x G M, / G i3(IR4.) and all integers n > 1, we have 



(2.5) E[L(x,/)]" = (2^)-"^^ j exp ^-i^Uj-x^Eexp ^i^UjX(tj)^ 



dudt. 



In the above and in the sequel, du = dui ■ ■ ■ dun and dt = dti ■ ■ ■ dtn- 
We start with the following moment estimates for the local time of X. 

Lemma 2.1. Let L = {Lt,t > 0} be the local time at zero of a strictly stable Levy process 
X = {Xt,t > 0} with values in M and index 1 < a < 2. Then for all < a < b < oo and 
all integers n > 1, 



.~ r(i-i)r(i + i + n(i-i)) \ 

^ fb-a\ A^l n\ {b - a)"(^-V°) 
- [~r) r(l + n(l-i)) ' 

where Ai > is the constant defined by ( f i. 7| ). Ln the case a = 0, we have the equality 
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Proof. Applying (|2.5|) with x = and / = (a, b] and making a change of variables Sj = tj — a 
(j = 1, • • • , n), we obtain 



E[|Lfe-Lfer] = (27r)- 

(2.8) 



n! 



(27r)' 

where 



" /" /" Eexp ( i VnjX(tj) i dud* 

/ / Eexp iy^UjX{sj + a) duds, 
Js„ Jm" \ / 



Sn = {{si, ■■■ ,Sn) ■■ < Si < S2 < ■ ■ ■ < Sn < b - a} . 

Let Vj = X^/Lj Ui (j = 1, • • • , n), then the last sum in (j2.8p can be written as 

n n 

(2.9) ^ -"j^lsj + a) = XI (^(^i + a) - ^(sj-i + a)) , 



where sq := — a so that sq + a = Denote 

m = lei 



tan(^) 



X 

Since the process X has stationary and independent increments, we have 



n 



n 



E exp i X^ Vj {X{sj + a) — X{sj^i + a)) 
(2.10) ^ 

" I "X^'^'^J' ~ 1 ( ~ isi + a)(j){vi)). 

It follows from (j2.8p . (|2.9p . (|2.10p and a change of variables that 



77 1 

(2.11) 'J (2 



X exp ( — (si + a)(p{vi)) dvds, 



where dv = dvi ■ ■ ■ dvn- 

We know from symmetry that for any s > 0: 



f f°° ( sv'^\ ( sv^u a-K \ 
/ exp ( — s (pM) dv = 2 exp ( I cos ( tan ( — ) I dv. 

Jm Jo \ X J \ X ' 2 V 



By a change of variable z = v", this equals 



2 1/^1 / s z\ ( szv .avTv , , 

«i„ »p (- -j CO. (- — .a„ {-))<;., 



which by Equation 3.944(6) in Gradsteyn and Ryzhik 22j equals 

,„,cos(larctan(.ta.(^))) . 

asi/<>(l + [^taii(ai)]2)''<'°> \a ^ 2 "J 

Combining the above with (|2.11l) we obtain 

„ n 

(2.12) E[\Lb - Lar] = n\ C7(a)" / {s, + a)"!/" His, - s,-_i)-i/" ds, 



where C(a) is the constant given by 

A/a-p(l 



x^/"r(i) /I 



(2-13) C{a) = + -ct- tan(f )) 

We denote the multiple integral in (|2.12p by J„. When a = 0, it can be evaluated in 
terms of the Gamma function. When a > the same induction method can still be applied. 
We include a proof for completeness. 

First we integrate over s„ € [sn-i, b — a] to get 



Next we integrate over Sn-i G [sn-2, & — a]. By changing variables twice v = Sn-i — Sn-2 
and t = v/(b — a — Sn-2), we obtain 

[Sn-l - Sn-2) ' " J OSn-l 

r / ^"'^"(^ - « - - Vf-^l'^dv 

1 - ^ Jo 

/■\(i-i/"M(l_t)(2-i/a)-i^i 

(6-a-.„,2)^(^-^/") r(i-^)r(2-i) 
1-^ r(i + 2(i-i)) 



(h-a-Sn-2?^^-^l^y 



r(i + 2(i-i)) 

Iterating this procedure, we derive 

Jn = TTT / (6-<!-«i)<"-'><'-"'»(si+o)-""<i»i 



pi. r(i + („-i){i-i))J„ 

^ T/i 1 \ri-l f(b-a)/b 

r(l + (n-l)(l-i)) Jo 

where the second equality follows from change of variables. 
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If a = 0, then the last integral equals Beta(l — ^, l + (n — 1)(1— ^)), where Beta denotes 
the Beta function. This and ()2.14p yield 

(2.15) Jn = , ^ " 1 X . 

r(l + n(l-i)) 

It follows from (f2l^ and (f2J5]) that for a = 0, 

(2.16) ^{\U-L^\\-^^^^-^^-^h^ IK 

By dLlI) and (f2l3]) . we have = r(l - C(a). Hence the desired result follows 
from (I2T6]) . 

If a > 0, then last integral in ()2.14p is an incomplete Beta function [22, 8.384(1) and 
9.100]. Write 

Ah-a)lh 

f{n, a,b)= / ^{n-m-l/a) (-^ _ ^yl/a 



(2.17) 

Then one can verify that 



f(n,a,b) <[-^) / t;("'^)(^-^/")(l-t;)"^/"dt; 

^^■^^^ '6 - "(i-i/")+i/"r(i - i)r(i + (n - i)(i - 1)) 



r(i + n(i-i)) 



and 

n(l-l/a)+l/a ^1 



(2-19) n(l-l/a)+l/a ^ 



& y l + (n-l)(l-i)' 

It follows from ([XT^ . (fmi) and that 

(-0) E[i..-..r] . (^)'" ^[;r;:;::!t;;" (^-<'-'^°'- 



Recalling = r(l - i) C7(a), we see that (j2.20p gives the upper bound in (j2.6l 
Similarly, we use (|2T^ . (|27ri]) and (|2T9]) to derive 



Nn(l-l/a) 



r22n EH^ / |ni-.f^-«V^° r(l-^)-^C(a)»n! 

(2.21) E[|L,-L.|]>^^j r(l + l + n(l-l)) " « 

which yields the lower bound in (j2.6p . This proves Lemma l2.1i □ 

Now we consider the moments of the increment Z^{h) — [a). As we will see in Remark 
14.11 the following lemma is sufficient for proving ()1.8p in Theorem 11.21 
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Lemma 2.2. Let = {W^{t),t G M} be a fractional Brownian motion of index H in 
M, and Lt be the local time at zero of a strictly stable process X = {Xt^t > 0} of index 
1 < a <2 independent of . Then for all < a < h < oo and all positive integers n, 



(2.22) Ci{n) {b - a)"(i-i/°) < E (\w"{Lb) - < C2{n) {b - a^^^-^ 



/a) 



where 



(2.23) Ci(n) = ^ / ( ^—^) (2^(2^) A,)- , "'^/^ + 

\/^ri-i V b J ^ ^ ri + i + n(l- 



V^r(i-i) \ b J ' r(i + i + n(i-i)) 

and 

In the above Ai > {) is the constant defined by ji. 7| ). Moreover, when a = we have the 
equality 

(2.25) E f |T^^(L,)r/^') = ^ (2^/(2^^) A^)" ^^'^(aTr + ^) ^nd-i/a)^ 

VI ^ ''^1 ; r(l + n(l-i)) 

Proof. Since Lt is a non-decreasing process, is iJ-self-similar with stationary incre- 
ments, and these two processes are independent, we use a conditioning argument to derive 

(2.26) E (|M^^(Lb) - W^^(L„)|"/^) = E (|W^^(l)r/^) E{\L, - LJ"). 



On the other hand, since W^{1) has a standard normal density, one can use a change of 
variables to verify that 

(2.27) E ( iW^iir/"") = 2"/(2^^) t(^ + -). 

\ / yTT \2H 2 J 

Combining I^Mi), with I^M) and ([22]) proves ([2:22]) and ([2:25]) . □ 

3. Analytic results: exponential integrability 

In this section we study the exponential integrability of the random variable (t) and 
some analytic properties of its logarithmic moment generating function. Our main results of 
this section are Theorems 13.11 and 13.41 which are the main ingredients for proving Theorems 
OandO 



Theorem 3.1. Let = {Z^{t),t > 0} be an a-stable local time H -fractional Brownian 
motion with values in M and 2H < a. Then for every 6 gM, 

(3.1) lim r logEexp (eZ^{t)) = Ai((9), 

2a 

where Ai is the function on R defined by Ai(^) = Bi Qa-^H fg^ all 6 (zW. Recall from lfL6\) 
and |j. that the constants Bi = Bi{H,a,x,'^) o-nd Ai are defined as 



(3.2) B^^"-^"' "^^ 



2a 
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and 

, _ r(l - ^)r(l)x^/"cos(^arctan(i.tan(^))) 
^ ' ' 7ra[l + (z.tan(^))2]i/(2") 

where v E [—1, 1] and X > are the constants defined in i2.1\) . 



Note that for 2H < a the above function Ai(-) is even, convex and difFerentiable on M, 
and the function A*(-) defined by (jl.Sp is the Fenchel-Legendre transform of Ai, that is, 
Al{x) = supegK {6x - Ai{6)) for all x £M. 

The proof of Theorem 13. II relies on explicit calculation of the moments of (t) and the 
following theorem in Valiron [40, page 44]. 

Lemma 3.2. Let f{z) = 'YlpLo'^p^^ entire function such that Cp ^ for infinitely 

many p's. For any r > 0, let M{r) = supj^j^^ |/(z)| . Then a necessary and sufficient 
condition for 

(3.4) hm i^^i^ = B 

r-^oo rP 

is that, for all values of e and all sufficiently large integers p, we have 

(3.5) —pcP/P<B + e, 

pe 

and there exists a sequence of integers pn, such that 

(3.6) hm ^ = 1, 

for which 

(3.7) lim —pnCp/P- =B. 

n—>oo pe 

Proof of Theorem \3.1[ Similar to the proof of Lemma 12.21 apply a conditioning argu- 
ment and the formula for the moment generating function of a Gaussian random variable 
to derive that for all € M, 

(3.8) Eexp(0Z^(t)) =Eexp (^^ t2H{i-i/a)^2H^ 

In order to prove (|3.ip . we show that, for 2H < a, the function f{z) = Eexp [zLf^) is an 
entire function and the coefficients of its Taylor expansion verify the conditions of Lemma 
[321 

Let us first consider the Taylor series 

(3.9) Mi(r) = y ^^^\ 

n=0 

We will make use of the following consequence of Jensen's inequality: For any constant 
7 > 1 and nonnegative random variable A, 

(3.10) (E(AW)y^^^^ < E(A^) < (e(aW+i;^''^^^^^^'^ 

Here and in the sequel, [7J denotes the largest integer < 7. 
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It follows from (|3.1U|) with A = Li and 7 = 2Hn, (|2.7|) in Lemma 12.11 and Stirling's 
formula that 



(3-11) ^1 - (^-^1 ^ 

where j4i is the constant in (jl.7p . In the above, x„ y„ means that, for all n large enough, 
Xn/Un is bounded from below and above by constant multiple of n~^. Here rj is a constant 
depending on H and a only. The omitted factors have no influence on the limit in (I3.12P 
below. 

By (|3.1ip . we see that the Taylor series in (|3.9p represents an analytic function on M if 
and only if 2H < a. In the latter case, we choose pi = ^"211 derive 

(3.12) lim ni^^ — '- ] =—(Ai ^ ' 



n^oo Pie \ nl J pi V (1 - ^)^"^/° 

Hence, Lemma 13.21 implies that 

(3.13) = P'^'"" y"'' 



r— ►oo 



rPi pi\ {I - i)i-i/a 



It follows from 1^ that Eexp {eZ^{t)) = Mi (f _ jjence 1^ follows from 

p.l3p and a simple change of variables. □ 

The proof of Theorem O shows that, if 2H > a, then E(e^^'^W) = 00 for all 9 > 0. 
Hence we cannot prove a large deviation principle for Z^{t) by applying the Gartner-Ellis 
Theorem. However, for studying the tail probability of (b) — Z^{a), it is sufficient to 
consider the exponential integrability of \Z^{b) — Z^{a)\^ for appropriately chosen f3 > 0. 

Proposition 3.3. Let < a < b < 00 be given constants. For any /3 > and t € M, let 

gp{t) = ^i^e*\w"{L,)-w"(La)\f'Y The following statements hold: 



(i) If < f3 < 2H%oi ' ^^sn the function gf3{t) is analytic on M. 

(ii) If [3 = 2H+a ' ^^^'^ dpi^) ^-^ analytic in (—00,(^0) for some 5q > 0. 

(iii) If P > 2H+a ' ^^^'^ dpi't) — ^ for all t > 0. 
Proof. Let us consider the Taylor series 

As in the proof of Lemma 12.21 we have 

E(|I^^(Lb) - VF^(La)|''") = E(|M^^(1)|'^") E(Lfe - L,|^^") 

^'•''^ =^2(^")/^r(^ + i)E(|L,-L.r'^). 

Applying (j3.10p to A = [L;, — La| and 7 = f)Hn and using (j3.14p . the moment estimates 
in Lemma |2. II and Stirling's formula, we derive 

E(r"(L.) - W"(K)r) ^ „(.-|-M) „-„(-f-M) 
n\ 
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where A2 is a constant which can be expressed exphcitly in terms of a, (3, H, Ai and b — a. 
It can be verified that ()3.15p imphes the conclusions in Proposition 13.31 □ 

It fohows from Proposition 13.31 that, for /? > 0, gp{z) {z E C) is an entire function if and 
only if /3 < 2H+a • ^'^^^ case, Theorem 13.41 further proves that gp{z) is of very regular 
growth in the sense of Valiron ( [iQ] ) . 

Theorem 3.4. Let be a fractional Brownian motion in M, and let Lt be the local time 
at zero of a strictly stable process Xt of index 1 < a < 2 independent ofW^. Then for all 
< P < ^^d. < a <b < 00, 

logEl exp (t\W" (Lh) - (La)\'^)] 

(3.16) lim ^ ^ ' — ^ ^' = ^3, 

where p = 2a-ap-2H(3 '^'^'^ ^3 ~ B-i{H,a,i',x, P) is the constant given by 

(3.17) B3 = i A^"' (b - <.)W->li-l/.) ( 0^"'"'Mmf" \ 



where Ai is the constant given in (1.1). 

Note that, when H = 1/2, we can take (3 = 1. In this case, is reduced to 

^ r(l - l)r(^)x^/"cos (larctan(..tan(n^))) " 
27ra[l + (i/tan(^))2]i/(2«) 

Proof of Theorem \3.4\ We start with the following elementary fact: Let H G (0, 1) be a 
constant. Then for all x > 0, 

In order to verify (j3.18p . first note that the first inequality holds for all x > because 
{[■^J,n > 0} is a subsequence of N, and the second inequality holds for all < x < 1. 
Hence it only remains to show the second inequality holds for all x > 1. This can be 
verified by grouping the terms in the expansion = 1 + Y.T=i TT ™ the blocks k e 
{L^^J +!,••• ) ^-iid noting that the number of integers in each block is at most 

1 + which is smaller than jj. 

Now let /3 e (0, 2^) be a constant and let t > 0. By taking x = t\W^{Lh) - W^{La)\'^ 
in (|3.18p . we have 

^ ^ti-mE{\w^{U)-w^{L,)fi-m) ^ 
1+ 2^ ^ 

n=l 

(3.19) , 



e 

< 



n=l 
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Denote the first and last terms in (|3.19|) by fi{t) and f2it), respectively. In order to prove 
(j3.16p . it suffices to show that for i = 1,2, 

(3.20) lim = B,, 

t-+oo tf^ 

where p = 2a-ap-2Hfi ^3 = B^{H, Q, V, X, (3) is given by p.l7p. 

This can be done by showing the coefficients of fi satisfy the conditions of Lemma 13. 2i 
Moreover, since the proofs for i = 1,2 are almost the same, we only prove p.20p for i = 1. 

Note that the coefficients Cp (p = 0, 1, • • • ) of fi{t) are given by 

_ E{\W^{Lb) - W^{La)\^^''/"i) 
if p = [-^J and Cp = otherwise. By Lemma 13.21 it suffices to show 

As in the proof of Lemma 12.21 we have 

By (j3.10p . the moment estimates in Lemma |2. II and Stirling's formula, we can verify that 
(3.23) H\W-iL,)-Wy.)fl-/-l) ^ ^LW^j ^.(i,_^-|) „-n(i,-^-|) 



where is a constant defined by 



(3.24) A3^<(^-ar"-^' ^ '^_,^,|-., 

Since {jj — — a) Hp = 1, we see that ()3.23p implies 

^ Lit J 1 = 



n—*oo 



where is given by p.l7p . This proves ()3.2ip and hence Theorem 13.41 □ 



4. Large deviations results: proofs of Theorems 11.11 and 11.21 

In this section, we first prove Theorems II. II and II. 21 Then we apply a maximal inequality 
due to Moricz, et al. [32] to derive upper bounds for the tail probabilities of the maxima 
maxtg[o,i] Z^{t) and maxtg[„^b] \Z^{t) - Z^{a)\. 
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Proof of Theorem \l.l\ Note that the function Ai(0) = BiO^-'^h in Theorem 13.11 is es- 
sentiaUy smooth and continuous on M. It follows from the Gartner-Ellis Theorem (cf. [201 

2H(q-1) 2H(Qf-l) 

Theorem 2.3.6]) that the pair (t "^-'^h Z^{t),t °'-^h j satisfies a large deviation principle 
with the good rate function 

K\{x) = sup {9x-Ki{e)), 

em. 

which is the Fenchel-Legendre transform of Ai. It is elementary to verify that AJ(j;) coin- 
cides with (|1.5p . This proves Theorem ll.il □ 



Proof of Theorem Let < /3 < 2H+a fixed. It follows from ()3.16p in Theorem [37 
and Davies' Theorem 1 in [T8| that 



an 1- ^oS^{\W''{L,)-W^iLa)f>u} -lu R^-i/(p-i) 

(^•1) ^MPT^ = -(l-P )ipB3) 

Here p = 2a-ap-2Hi3 ■ Letting x = u^/^ and simplifying the right-hand side of (14. ip . we 
obtain 



logP{|VF^(Lfe)-W^^(L,)| >x} 
lim 



(4.2) 




a + 2H HA? \ „ 



2a 

This finishes the proof of Theorem 11.21 □ 

Remark 4.1. Theorem 11.21 can also be proved by using Lemma 2.3 in Konig and Morters 
[29j (note that their assumption p G N can be replaced by p > 0) and the moment results in 
Lemma 12.21 The proof of Lemma 2.3 in [29] is based on a change-of- measure technique in 
large deviations. We remark that Lemma 2.3 in [29] is equivalent to Corollary 2 in Davies 
|18j . hence it can also be proved by using an analytic method. 



Similar to the proof of Theorems 13 . 1 1 and 1 1 . H we obtain the following theorem for the local 
time Lt at zero of X. They are in complement to the results of Hawkes |23j and Lacey |30| 
on the tail asymptotics for the local time Lt and the maximum local time max^^gK L(x, t), 
respectively. 

Theorem 4.2. Let L = {Lt,t > 0} be the local time at zero of a real-valued strictly stable 
Levy process X = {Xt^t > 0} of index 1 < a < 2. Then the following two statements hold: 
(i) For all < a < b < oo, 

/A o\ r logIE[exp(t {Lb - Lg))] ^ x 

(4.3) lim = [b — a) C{a,iJ,x), 

where C{a,i',x) is the constant defined by 

>(1 - i)r(i)x^/"cos (iarctan(z.tan(^)))" " 



C{a,i/, x) 



7ra[l + (z.tan(^))2]^/2" 
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(ii) The pair {t '^^^^ ^^(L;, — L^), i"/^" satisfies LDP with good rate function K2{x) = 

-(a-l) 

if X > and A2(x) = cxd if x < 0. That is, for every 



{b-a)C{a,iy,x} 
Borel set F C M, 



(4.4) 

and 
(4.5) 



limsupt logPjf (Lb - -La) G f| < - in£A2(a;) 



liminf t «-i logP<^ t 



>oo 



1 

" Q-1 



(Lb - L,) G F > 



} 



inf AS(x). 



Proof. Similar to the proof of Theorem 13. H Eq. (j4.3p fohows from Lemma 
moment estimates in Lemma |2.1[ 
It fohows from (i) that for ah > 0, 

logE[exp(0t(Lfe-L,))] 



and the 



(4.6) 
Denote 
(4.7) 



hm 

4— »oo 



(6-a)C(a,i.,x)e"/("-'^ 



{b-a)C{a,u,x)0''/^''-^^ if ^ > 0, 

if < 0. 



Then, in the terminology of [20], A2 is an essentially smooth, continuous function on M and 
its Fenchel-Legendre transform is given by 

^)(6-a)C(a,z.,x) 



(4.8) 



AUx) 



00 



if X > 0, 
if X < 0. 



Therefore, as in the proof of Theorem ll.il Part (ii) follows from (14.61) and the Gartner-Ellis 
Theorem. □ 

Remark 4.3. Let F = [l,oo) in (14. 4p and (14. Sh . we obtain the tail probability 



(4.9) 



lim 

a;— >oo 



logP{Lfe-L„>x} 



1 

a 



a 



a 



1 



{b-a)C{a,u, x) 



-(a-l) 



In case a = 0, the limit in (|4.9p is weaker than the best known. By using the connection 
between Lt and a stable subordinator of index 1 — 1/a, one can derive more precise result 
(i.e., without the logarithm) on limiting behavior of PjL;, > x} as x ^ 00. See Hawkes 
[23]. 

Theorem 11.21 is concerned with asymptotic behavior of the tail probability P{|Z^(6) 
—Z^{a)\ > x} as X ^ CXD. In many applications, however, it is useful to have sharp bounds 
on¥{\Z^{b)-Z^{a)\ > x} and¥{maXa<t<b\Z^ (t) - {a)\ > x} that hold for all x > 0. 
In the rest of this section, we consider these questions and in the next section we use them 
to derive upper bounds for the local and uniform moduli of continuity for Z^ . 

Lemma 4.4. There exists a finite constant > 0, depending on H, a, v and x only, such 
that for all < a < b < 00 and all x > 0, 

^2a/{a+2H) \ 



(4.10) 



'{iZ^ib) - Z^{a)\ > x} < exp -A4 



(5 _ a)^H{a-l)/{a+2H) 
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Proof. We consider the random variable 

\Z^{b) - Z"{a) 



A 



(5 - a)^("-i)/" ■ 

As in the proof of Lemma 12.21 we apply a conditioning argument and Lemma 12.11 to show 
that for all integers n > 1, 

(4.11) E(A") < Al n"("+2^)/(2°) , 

where > is a constant depending on H^a,v and x only. 

For any constant > 0, the Markov inequality and (14. lip imply that for all u > 0, 

An n{a+2H)/(2a) 

r4i2^ ° 



By taking the constant Aq > eA^ and n = [u^^^ \ , we obtain 
(4.13) F{\Z^{b) - Z^{a)\ > (6 - a)-f^(°-^)/" u} < exp /"--^^/C^+s^^) ' 

It is clear that (fiJO|) follows from (fiJ3|) by letting x = (6 - a)'^("~^)/'*'u. □ 

Next we apply Lemma 14.41 and a result of Moricz, et al. [32j to prove the following 
theorem. 

Theorem 4.5. There exist positive and finite constants A-j and Ag, depending on H,a,v 
and X only, such that for all < a < b < 00 and all x > 0, 



^2a/{a+2H) 

-As — 



)2H{a-l)/{a+2H) 



(4.14) p| max \Z^{t) - Z"{a)\ > x 1 < yly exp ( - 

a<t<b J \ 

Proof. For any integer n > 2, we divide the interval [a, b] into n subintervals of length 
(b — a)/n. Let tn,i = a + ^^^^ [i G {0, 1, . . . ,n}) be the end-points of these subintervals. 
By the sample path continuity of Z^ , it suffices to show that 



(4. 



15) p| max \Z^{tn i) - Z^ {a)\ > x 1 < A7 exp ( -Ag 

l<i<n ' J I 



X 



2a/{a+2H) 



iyH{a-l)/{a+2H) 



for all integers n > 2. 

To this end, we define the random variables = Z^ {tn,i+i)— Z^ {tn,i) for i G {0, 1, 
1}. Then for all integers < j < < n, we have 



, n- 



k-l 



(4.16) ^^(tn,fc) - ^^(tn,i) =Y.^i--= S{j, k). 



Applying Lemma 14.41 we see that for all integers j < k, 



(4.17) r{\Sij,k)\ >x}< exp ^-^4 ( 



2H{a-l)/ia+2H) 
!^ \ 2a/{a+2H) 

{h-a)[k-j)' 
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Using the notation in [J2J, we denote (f){x) = x^"/^"^'^^) and 

/(h-n)(k- i) \ ■iH{a-l)/(a+2H) 

For simplicity denote r = 2H [a — l) / {a + 2H) . Since r G (0, 1), the concavity of the function 
1 1— > f imphes that for ah integers l<i<j<k<n, 

(4.18) gii,j)+g{j + l,k)<2^-''g{i,k). 

Hence the function g satisfies the property of quasi- superadditivity with index Q = 2^"*" 
in [32]. Moreover, the functions (p and g satisfy ah the other conditions of Theorem 2.2 
in Moricz, et al. [32]. Consequently, the latter implies the existence of positive and finite 
constants and ^8 (depending on H, a, v and x only) such that (|4.15|) holds. This proves 
Theorem 14.51 □ 

Remark 4.6. Note that, when H = 1/2, one can apply the refiection principle of Brownian 
motion and conditioning to prove Theorem 14.51 Our method is much more general. 

5. Applications 

Applying the large deviation results in the previous section, we establish uniform and 
local moduli of continuity for . 

Theorem 5.1. Let Z" = {Z^{t),t > 0} be an a-stahle local time H -fractional Brownian 
motion with values in R. Then there exists a finite constant ^9 > such that for all 
constants < a < b < 00, we have 

(5.1) limsup sup sup ■ -, — „ ' , < ^dg a.s. 

/^iO a<t<r/xO<A/iH(a-l)/a(logl//,)(°+2^^)/(2") - 

Proof. For every t > and /i > 0, it follows from (14.131) that 

(5.2) ¥{\Z^{t + h) - Z^{h)\ > Aeh^'^'^-^^/'^u} < exp (-n2"/(°+2^) 

Hence Z^ = {Z^{t),t > 0} satisfies the conditions of Lemmas 2.1 and 2.2 in [I3| with 
cr{h) = and /3 = ^^2H ' Consequently ()5.ip follows directly from Theorem 3.1 in 

[13]. " □ 



Csaki, Foldes and Revesz [T7] obtained a Strassen type law of the iterated logarithm 
(LIL) for Z(t) = W{Lt) when Lt is the local time at zero of a symmetric stable Levy 
process (see Theorem 2.4 in [17]). Part (i) of the following theorem extends partially their 
result to Z^ and Part (ii) describes the local oscillation of Z^ in the neighborhood of any 
fixed point. 

Theorem 5.2. Let Z^ = {Z"{t),t > 0} be an a-stable local time H -fractional Brownian 
motion with values in M. The following statements hold: 
(i) Almost surely, 

maXo<^<t |^^ (s)| ^ ^-(a+2H)/{2a) 

t^oo^ i-f^(o^-l)/°(loglogt)* 
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(5.3) lim sup — — — — (,+2f/)/(2») ^ ^8 



(ii) For every t > 0, almost surely 

^ ( log log 1/h) - ' 

In the above, is the constant in 

Proof. Since both (j5.3p and (|5.4p follow from Theorem 14.51 and a standard Borel-Cantelli 
argument, we only prove (|5.3p . 

Fix two arbitrary constants 7 > A^^ and p > 1. For every integer n > 1, let = p"" 
and consider the event 

En=L: max > (-D/" [/(Tjj , 

where C/(t) = (7 loglogt)(°+2H)/(2a)^ j^. fono^g f^om Theorem US] that 

(5.5) P(ii;„) < exp i^-As ^ " j < ^10 n-^«^ 

Since Agj > 1, we have Yl'n'=i^i'En) < 00. The Borel-Cantelli lemma implies that 

(5.6) limsup^^^^;?^ffL^^<l a.s. 

Note that Tn+i/Tn = p for every n > 1. We use the monotonicity to derive that for all 
t G [T„, r„+i], 

, . maxo<^<t |^^(s)| maxo<s<T„+i {s)\ U{Tn+i) 



Equations (j5.6p and (|5.7p imply 

f 5 8) lim SUD maXo<,<f |Z-^(s)| ff(a-l)/a (a+2H)/{2a) 

(^O.OJ nmsup ^ (a+2H)/(2a) - ^ ^ 

We obtain (j5.3p from (j5.8p by letting 7 | Ag ^ and p j 1 along rational numbers. □ 

Remark 5.3. We believe that, up to a constant factor, both uniform and local moduli of 
continuity of Z^{t) are sharp. However, we have not been able to prove this due to the 
lack of information on the dependence structure of the process Z^ . 
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